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ABSTRACT 

A large body of evidence has demonstrated that the global rest-frame optical and in- 
frared colours of galaxies correlate well with each other (i.e. {u — g), with (r — 8/im)), 
as well as with other galactic properties such as surface brightness and morphology. 
However the processes that lead to the observed correlations are contrary; the stellar 
light that contributes to the optical is readily absorbed by dust which emits in the IR. 
Thus on small scales we expect these correlations to break down. We examine here 
seven nearby galaxies ranging from early- to late-types, which have all been smoothed 
to the same physical scale and signal-to-noise ratio, using data from the optical to 
the mid-IR (ix-band to 8/im). By examining these galaxies on a pixel- by-pixel basis 
we demonstrate that there is disconnect between the optical and IR when normalized 
to the near-IR (i^-band). For five of the seven galaxies we can decompose this dis- 
connect into two distinct components through a Principal Component Analysis of the 
i7-band normalized spectral energy distribution of the pixels: one mainly correlated 
with variations in the IR, the other correlated with variations in the optical. The two 
exceptions are the elliptical galaxy NGC4552, whose SED can be well reproduced us- 
ing a single principal component, and the highly inclined spiral NGC3521, due to its 
complex dust geometry. By mapping these two components in the five "regular" galax- 
ies, it is clear they arise from distinct spatial regions. By comparing these components 
with the surface brightnesses in Ha and i^-band we demonstrate that the IR dom- 
inated component is strongly associated with the specific star-formation rate, while 
the optical-dominated component is broadly associated with the stellar mass density. 
However, when the pixels of all galaxies are compared, the well known optical-IR 
colour correlations return, demonstrating that the variance observed within galaxies 
is around a mean which follows the well-known trend. As a final step, we extend this 
work by examining the extremely tight correlations observed between the IRAC-near- 
IR colours, and demonstrate that these correlations are tight enough to use a single 
IRAC-near-IR colour (i.e. Sjj.m — H) to determine the fiuxes in the other IRAC bands. 
These correlations arise from the differing contribution of stellar light and dust to the 
IRAC bands, enabling us to determine pure "stellar" colours for these bands, but still 
demonstrating the need for dust (or stellar) corrections in these bands when being 
used as stellar (dust) tracers. 

Key words: galaxies: general, photometry, stellar content, ISM; infrared: galax- 
ies; galaxies: individual: NGC3521, NGC4254, NGC4321, NGC4450, NGC4536, 
NGC4552, NGC4579. 



1 INTRODUCTION 

One of the remarkable aspects of astronomy is the strong 
correlations that exist between the broad observable prop- 
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erties of galaxies. The morphological sequence of galaxies 
first noted by [Hubble (1926) correlates well with surface 
brightnesses, sizes, and optical colours. These in turn cor- 
relate with and arise from the intrinsic physical properties 
of these galaxies, such as the total stellar masses (M*), star 
formation rates (SFR), and mean stellar ages and metallici- 
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ties (for reviews on these correlations see Roberts & Haynes 
T994| and |Kennicutt|[T998t 



The rest-frame optical colours of galaxies broadly cor- 
relate with each other (i.e. (u — g), with (g — r)), as well as 
with surface brightness. While forming continuous sequence 
in colours, galaxies can be broadly classified into 'blue' and 
'red' galaxies, especially when compared with luminosity or 
absolute magnitude, with a strong correlation with late- and 
early-type respectively (Blanton et al. 2003 ). The variation 
in optical colours arise from differences in the star forma- 
tion histories of galaxies, with UV to blue light dominated 
by short-lived massive stars, while the red to near infrared 
light is more sensitive to the total amount of stars. This 
is complicated however by the presence of interstellar dust. 



which attenuates and reddens the stellar light (see e.g. Ken- 
nicutt|1998||Bruzual Char lo t |2003t . Thus the colours of a 



galaxy broadly measure the ratio of young to old stars, which 
is generally parameterized by the ratio of current to past star 
formation (b =SFR/ <SFR>), or the ratio of the current star 
formation rate to the total stellar mass, the specific star for- 
mation rate (sSFR=SFR/M*). In the local universe, these 
intrinsic physical parameters, and thus a galaxy's colours, 
are related to the total stellar mass of a galaxy (e.g. Brinch- 



mann et al.|2004 ), with more massive galaxies having lower 



sSFRs and thus redder colours. The stellar mass, and thus 
the sSFR and colours, are related to the size, concentration, 
surface brightness, and mean metallicity of a galaxy (e.g. 



Kauffmann et al.|[2003l [Tremonti et al.|2004| [Gallazzi et aT] 



2005[ 



In a similar manner, the IR colours of galaxies are also 
strongly associated with the intrinsic properties of the galax- 
ies, with the IR fluxes increasing relative to the near-IR 
fluxes (i.e. becoming redder) from early- to late-type galax- 
ies (see e.g. Kennicutt||1998 especially fig. 4), and with in- 
creasing sSFR ( |da Cunha et al.|[2008 ). As with the optical 
colours, the infrared is sensitive to the younger stellar pop- 
ulations, as the IR emitting dust is heated predominantly 
by UV-blue light. In fact, the fluxes in both of these wave- 
length regions can be used as measures of a galaxy's SFR 
(Kennicutt 1998) and a clear correlation between the global 



optical and IR colours is observed, in the sense that optically 
blue galaxies are also found to have enhanced IR emission 



(e.g. |Hogg et al.|2005 ). 

However on a physical basis, the processes that lead 
to the fluxes in these two wavelength regimes are contrary: 
the radiation from young stars that dominate the UV-blue 
optical colours is also that preferentially absorbed by dust 
which then re-emits in the IR (see e.g. [Dra ine| |2003 l . This 
contrary nature is exacerbated by the association of the 
youngest stars with the clouds of gas and dust from which 
they form, leading to relatively higher attenuation observed 



in star forming regions (see e.g. Calzetti 



1997 



Chariot & 



[Fall 2000 ) . This contradiction of the observed galaxy scale 
correlation of both blue colours and IR excess and the phys- 
ical contrary nature of dust absorption and UV-blue light 
emission indicate that on some spatial scale within galaxies 
this correlation must break down. 

In some respects this difference is expected as the blue 
light and infrared emission arise from physically distinct 
components (i.e. stars and dust), and this correlation will 
break down once these different components are resolved 
out, as seen within our own Galaxy. However, this distinc- 



tion between IR and blue light exists on larger scales within 
galaxies, as we show here by examining seven nearby galax- 
ies that span a range of galaxy types from early to late. 
By examining all galaxies on the same spatial scales we 
demonstrate that the variation of the IR light and the blue 
light occurs in spatially distinct regions, and show the scales 
over which these regions occur. Importantly, these distinct 
components appear to be associated with different physical 
quantities that are correlated on galaxy scales: the specific 
star formation rate (sSFR) and the stellar mass density. 

In Section [2] we introduce the sample of galaxies, the 
mult i wavelength dataset and the relative reduction. In sec- 
tion [3] we analyse the resolved SEDs of the individual galax- 
ies, showing the correlations that exist to different degrees 
between the IR and optical colours, analyzing them by 
means of principal component analysis, and showing how 
they depend upon other local observables, namely the NIR 
and Ha surface brightness. In Section [4] we show how these 
local colour correlations (and lack thereof) result in stronger 
global correlations. Section [5] is devoted to an in-depth in- 
vestigation of the SEDs between 1.65 and 8/im, which we 
describe as a universal 1-parameter family: fitting formulae 
to obtain the luminosity at intermediate wavelengths given 
H-hand and S/nm and estimates of purely stellar SEDs are 
provided. Section [6] concludes the paper with a summary of 
the principal results. 



2 SAMPLE AND IMAGING DATA 

The present sample of 7 galaxies is drawn from the origi- 



nal sample of nine galaxies analyzed in Zibetti et al.| ( 2009 



ZCR09 hereafter). We require galaxies to be part of the 
Spitzer Infrared Nearby Galaxies Survey (SINGS, Kenni- 
cutt et al. 2003 ) , a comprehensive imaging and spectroscopic 
study of 75 nearby galaxies (D < 30 Mpc) conducted in the 
IR with the Spitzer Space Telescope, and have SDSS imag- 



ing (York et al. 2000) as well as deep NIR imaging in H band 



(1.65/im) either from GOLDMine (Gavazzi et al. 2003) or 
UKIDSS ( Lawrence et al.||20 07). In addition" we require the 
galaxies to be at most ^17 Mpc away (corresponding to 
the assumed distance to the Virgo cluster, 'Gava zzi et al.| 
1999), in order to resolve scales of approximately 200 pc 
at all wavelengths, from SDSS ix-band to IRAC 8/im. This 
leaves us with a small, yet representative, sample of regular 
galaxies in the local Universe. The 7 galaxies are listed in Ta- 
ble[l]along with their coordinates (according to NED/RC3), 
assumed distances, corresponding angular scale, source and 
depth of the NIR imaging (see below). With these seven 
galaxies we cover the full morphological range of regular 
"giant" (3.9 10^ < M^/Mq < 9.3 10^°) galaxies, from ellip- 
tical to Sc spiral, also corresponding to a range of specific 
star formation rates from < 10~^^ to ~ 10~^ yr~^. In Table 
|2] we present these and other essential morphological and 
physical properties of the galaxies in the sample. The in- 
ferred metallicities also span the full range for galaxies of 
such stellar mass. Only two spirals have an apparent incli- 
nation in excess of 45 degrees, thus allowing us to study 
not only a "clean" sample of low inclination discs, for which 
the radiative transfer in the dusty ISM is expected to be 
relatively simple, but also to explore possible complications 
due to a larger inclination. We note that nuclear activity 
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Table 1. The sample. 



Name 


RA 


Dec 


Distance 


Angular scale 


NIR source 


Limiting SB{H) 




(J2000.0) 


(J2000.0) 


Mpc 


pc arcsec"^ 




MJy Sr-i 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC3521 


Ilh05m48.6s 


-00d02m09s 


9.2 


45 


UKIDSS 


0.30 


NGC4254 


12hl8m49.6s 


+14d24m59s 


17.1 


82 


GOLDMine 


0.21 


NGC4321 


12h22m54.9s 


+15d49m21s 


17.1 


82 


GOLDMine 


0.32 


NGC4450 


12h28m29.6s 


+17d05m06s 


17.1 


82 


GOLDMine 


0.35 


NGG4536 


12h34m27.0s 


+02dllml7s 


17.1 


82 


GOLDMine 


0.39 


NGC4552 


12h35m39.8s 


+12d33m23s 


17.1 


82 


UKIDSS 


0.47 


NGC4579 


12h37m43.5s 


+Ild49m05s 


17.1 


82 


GOLDMine 


0.53 



from AGN does not influence our analysis significantly, as 
we include only weak or obscured AGNs (no type 1) which 
affect only a very limited number of central pixels. The sta- 
tistical effect of these AGN affected pixels over the several 
thousands of pixels in each galaxy is negligible. 

The SDSS and if-band images are processed, mosaiced 
(when required) and calibrated as explained in ZCR09. 
Spitzer IRAC images are publicly available from SINGS, 
already calibrated in MJy/Sr and sky-subtracted. For part 
of our analysis we also use narrow-band Ha imaging, pro- 
vided as an ancillary dataset of SINGS: details about the 
flux calibration of the "net" (i.e. continuum-subtracted) im- 
ages is given in the technical documentation that accom- 
panies the data release. No Ha image is provided for the 
elliptical galaxy NGC4552, which allegedly does not have 
any detectable HII regions. 

For all images the properties of the background, which 
is the dominant contributor to the noise budget at all wave- 
lengths, are carefully computed. The mean background is 
subtracted by fitting a linear polynomial to the regions sur- 
rounding the galaxies, which are selected not to display any 
strong contamination from bright sources nor to be affected 
by obvious artifacts. Noise properties in terms of pixel-based 
r.m.s. and fluctuations on scales of the order of a fraction 
of the galaxy's effective radius are also determined in these 
regions and are then used to determine the adaptive smooth- 
ing strategy and the surface brightness cuts to be applied in 
the following analysis. 

We manually remove all obvious contaminating sources 
(foreground stars, background galaxies) by replacing those 
pixels with the surrounding average background plus noise, 
using the task imedit in IRAF. For more extended or una- 
mendable artifacts (e.g. the extended blooming due to the 
bright nucleus of NGC 4536 in IRAC bands) we simply mask 
out those pixels (and the corresponding ones in all bands) 
and discard them from the following analysis. 

Pixel- by- pixel SED analysis requires two more steps: 
the placement of all images on a common plate scale, and 
to ensure that in each pixel at all wavelengths a sufficiently 
high signal-to-noise ratio (S/N) is provided. The imaging 
data involved in this analysis have point spread functions 
(PSF, for optical or if-band) or point response functions (for 
Spitzer IRAC bands) FWHM ranging from approximately 
1 arcsec (some SDSS images in good seeing conditions) to 2 
arcsec (IRAC S/xm), while pixel scales range from 0.4 (SDSS, 
UKIDSS) to 1.6 (some GOLDMine if-band images) arcsec 
per pixel. By re-sampling the images with a common pixel 
scale of > 2.5 arcsec pix~^ we ensure that most of the flux of 



a point-like source is contained in one pixel. We can there- 
fore avoid convolving all images to a common PSF, which 
would be hardly feasible for a number of H-band images 
where the PSF is severely under-sampled and not enough 
stars are available to compute it properly. We choose there- 
fore to resample the images of the 6 galaxies belonging to the 
Virgo cluster to a pixel scale of 2.5 arcsec, corresponding to 
207 pc pix~^ for the assumed distance to Virgo of 17.1 Mpc 
( Gavazzi et al.|1999 ). For the remaining galaxy, NGC 3521, 
we resample to a pixel scale of 4.65 arcsec pix~^ , correspond- 
ing to the same physical resolution of 207 pc pix~^ as for 
the other galaxies. 

In order to extend our analysis with a sufficiently high 
signal-to- noise ratio (S/N) out to a large fraction of the op- 
tical radius, as in ZCR09, we perform an adaptive median 
smoothing of the images using the code ADAPTSMOOTH 
fZibetti 2009 ) , which also allows the matching of the variable 



smoothing kernel among all bands. A minimum S/N=20 per 
pixel is required in all bands, except ix-SDSS, for which a 
minimum S/N of 10 is deemed as sufficient. In fact, this band 
is imaged with a significantly lower sensitivity in the SDSS, 
such that a more massive smoothing would be required and 
a smaller radial extent could be covered if the same S/N=20 
had to be required. On the other hand the large leverage of 
this band on the SED makes it possible to study pixel-to- 
pixel variations even with an accuracy of 10% only, as we 
show in the following analysis. We note that the Ha images 
do not contribute to determine the size of the smoothing ker- 
nel required at each pixel (as this would imply missing a sig- 
nificant fraction of pixels where no substantial Ka emission 
is observed), but they are nevertheless smoothed according 
to the common kernel as for all other bands. Background- 
dominated noise is assumed for all images, with a uniform 
r.m.s. equal to the one measured in the background regions. 
A maximum smoothing radius of 4 pixels (2 for NGC 3521) 
corresponding to ^ 800 pc is allowed. All pixels for which 
the surface brightness (SB) cannot be computed with the 
minimum required S/N even after smoothing over the max- 
imum radius are flagged and discarded from the rest of the 
analysis. Furthermore, we also discard pixels whose SB is 
less than 10 times the typical large scale background fluctu- 
ations, in order to avoid spurious effects at the low SB end 
of the pixel distribution. The limiting SB (if) that is reached 
as a result of these cuts is given in column 8 of Table [l] 

Images of the galaxies in the ten broad bands and in Ha 
after editing, adaptive smoothing and SB cuts are shown 
in Fig. |A1| to |A7| in Appendix ^ coded in grey intensity 
according to SB in units of MJy Sr~^. In the rest of the 
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Table 2. Physical properties of the sample galaxies. 



Denomination Morph. type Inclination M* SFR sSFR Z(gas) Nuclear classification 

degrees log M© log M© yr-^ log yr-^ 12 + log (O/H) 
(1) (2) (3) (4) (5) (6) (7) (8) 



NGC3521 


SAB(rs)bc 


64 


10.52 


0.33 


-10.19 


9.01 


LINER 


NGC4254 


SA(s)c 


30 


9.71 


0.82 


-8.89 


9.17 


LINER 


NGC4321 


SAB(s)bc 


32 


10.46 


0.55 


-9.91 


9.11 


LINER 


NGC 4450 


SA(s)ab 


43 


10.80 


0.03 


-10.77 


9.13 


LINER 


NGG4536 


SAB(rs)bc 


67 


9.59 


0.49 


-9.10 


9.00 


NonJ«] 


NGC 4552 


EO-1 




10.97 






9.12 


LINER 


NGC 4579 


SAB(rs)b 


38 


10.06 


-0.06 


-10.12 


9.22 


LINER 



Notes. Columns: (1) NGC name. (2) Morphological type according to NED/RC3 ( |de Vaucouleurs et al.liggl] ), following the classification 
scheme of |de Vaucouleurs] ( | 195 9 1 ): S indicates spiral, E elliptical followed by a number indicating the ellipticity (0 means perfectly round); 
A is normal unbarred spiral, B barred, AB weakly barred; (s) is s-shaped, (r) indicates the presence of a ring, (rs) both ringed and s- 
shaped; a, b, c, d are the same as in the classical Hubble spiral classification. (3) Inclination (for spirals only), according to de Vaucouleurs] 
]et al.| ( ll99l] >; |Jarrett et al.| ( l2003l >. (4) Stellar mass from Skibba et al. (ApJ in press) assuming a |Kroupa| ( ]200lt IMF, except NGC 4450 
and NGC 4552. Stellar masses for these two galaxies have been computed multiplying the total if -band luminosity given by 2MASS-LGA 
( jJarrett et al.]2003| ti mes the M* /Ljj computed according to ZCR09 using the g — i, i — H colours measured from our images and 



corrected to a 



Kroupa| ( l200ll> IMF. (5) Sta r formation rate from Skibba et al. (ApJ in press), based on Ho; + 24/im luminosity, except 



for NGC 4450 which is from |Calzetti et al.](l2010t. (6) Sp ecific star formation rate, sSFR=SFR/M*. (7) Gas phase metalhcity, from the 
characteristic oxygen abundance from [Moustakas et al.] (|2010| based on the [Kobulnicky Kewley] ( l2004| calibration, except for NGC 
4450, NGC 4552, and NGC 4579, which are use the luminosity-metalhcity relation. (8) Nuclear activity, from NED. 



^ Moustakas et al. ( 2010) > find line ratios consistent with very low-level AGN activity. 



paper we discuss the analysis and the results for the full 
sample, but we show specific figures and plots only for two 
galaxies which are representative of late star- forming spirals 
(NGC 4254) and early type galaxies (NGC 4450), with the 
remaining five being presented in online appendices. 



3 SPATIALLY RESOLVED SEDS: INTERNAL 
TRENDS WITHIN GALAXIES 

3.1 Pixel-by-pixel SEDs and colour-colour 
relations 

The processing described in the previous Section provides 
us with an optical to mid-IR SED for each valid pixel, de- 
scribed by an array of values for the SB in MJy Sr~^ in 
ten broad bands (u^ g, r, i, z, H, 3.6, 4.5, 5.7 and 8/im) 
plus the surface brightness in erg cm~^ sec~^ arcsec"^ for 
the Ka line emission. At each pixel we normalize log lyfu in 
each band to H band in order to be able to compare on the 
same scale the SED of pixels spanning about two orders of 
magnitude in SB and hence to study the spatial variations 
of the SED shape across the galaxies. 

The choice of H band is dictated by three main rea- 
sons which make it the natural point of normalization: i) 
1.65/im divides the SED into one part of pure stellar emis- 
sion (possibly dust extincted), shortward of it, and one part 
affected by dust emission, longward of it. ii) the H band is 
the one that minimizes the effects of dust both in emission 
(see e.g. Fig. 5 of |da Cunha et al.|2008 ) and in absorption, as 
the extinction is always smaller at longer wavelengths. The 
z band is equally unaffected by dust emission, but is more 
sensitive to extinction; on the other hand, 3.6/im, although 
even less affected by extinction than the H band, is severely 
affected by dust emission, as we show below, in) if-band 
luminosity is probably the best proxy to stellar mass, as 



shown e.g. in ZCR09, showing the minimum dependence of 
M/L on age, metallicity and dust; therefore normalizing to 
H-hand is close to normalizing to the stellar mass density at 
each pixel, a much more fundamental parameter than any 
surface brightness. 

The distribution of pixels in SED for each galaxy is 
given for reference in the bottom right panels of Figs. |A1| 
to |A7| respectively, in the online Appendix |A] The num- 
ber of pixels per bin of normalized logiyfu in each band is 
represented by the colour coding. Galaxy SED distributions 
display different amount of scatter, which vary upon mor- 
phology and wavelength considered. As we show in more 
detail in the following analysis, most of the scatter is con- 
tributed by the mid-IR and the u bands. It is interesting 
to note that the integrated SED of a galaxy (flux-weighted 
average of pixels, signs in the figures) does not perfectly 
correspond to the (number- weighted) mean pixel SED ('x') 
or the most frequent SED (peaks of the distributions). A 
similar effect was already noted by ZCR09, where it turned 
out to cause a significant bias in stellar mass estimates, and 
can be simply explained by the fact that regions of different 
brightness are often characterized by systematically differ- 
ent sed£] 

As a first step, we study the correlations between the 
normalized fluxes in all possible pairs of bands for all pix- 
els within each galaxy. This is shown in Fig. ^ and [2] for 
NGC 4254 and NGC 4450, respectively, with the remain- 
ing galaxies shown in the online figures |B1| to |B4[ In each 
panel of the figures we plot the distribution of pixels in 



^ A typical example occurs in spiral galaxies with a bright central 
bulge: the pixels dominated by the red bulge light are fewer but 
brighter than the more numerous but dimmer pixels in the disc 
where the SED is bluer. As a result, the flux-weighted SED will 
be redder than the number-weighted one. 
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Figure 1. Distribution of the pixels of NGC4254 in the "colour-colour" spaces defined by the logufjy fluxes normalized to if-band, for 
the u, g, r, i, z, 3.6, 4.5, 5.7 and 8/im bands. The grey-scale intensity linearly traces the density of pixels at each point in the colour-colour 
spaces. The Spearman rank correlation coeflicient for the distribution of pixels is written at the bottom right corner of each corresponding 
panel. The histograms on the diagonal show the distribution of pixels in normalized log ufi, colour space as labeled on the y-axis. The 
total number of pixels analyzed in NGC 4254 is given at the top of the flgure. 



the (logUYfuY - ^OgUHfiyn) ^S. (logUx fi^x - ^OgVHfvn) 

space (where X and Y are any two broad bands among 
g, r, i, z, 3.6, 4.5, 5.7 and S/xm), which is effectively 
a colour-colour space {Y — H) vs. {X — H), modulo a 
linear transformation. The panels on the diagonal display 
the distribution of the fraction of pixels as a function of 
{log uyfuY ~ logUH fun)^ which allows to directly derive a 
more quantitative information about the relative contribu- 



tions of pixels at different colour. In all following discus- 
sion we refer to ([oguxfux — ^ogUH fun) "[-^] colour" for 
brevity. 

By comparing Fig. ^and|2]it is immediately clear that 
the late-type, star-forming galaxy NGC 4254 spans a much 
broader range of colours than the early type NGC 4450, 
which is overall very uniform. Apart from this, both galaxies 
display common features in the colour-colour distributions. 
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Fraction 
0.0 0.15 0.30 



4692 pixels 



Fraction 
0.15 



Fractio 
0.0 0.1 



1.2 -0.8 -0.4 0.0 -0.6 -0.3 0.0 0.3 -0.25 0.00 0.25-0.15 0.00 0.15 0.30 0.00 0.15 0.30 -0.8 -0.6 -0.4 -1.2 -0.9 -0.6 -1.5 -1.0 -0.5 0.0 
[u] [g] [r] [i] [z] [3.6/.nn] [4.5/.nn] [5.7^nn] 



Figure 2. Same as Fig. [T] but for NGC4450. 



such that one can divide up these figures into three distinct 
regions, corresponding to different spectral correlations. In 
the upper left panels we observe very significant correlations 
between optical colours, from [u] to [z]. Similarly, the bottom 
right panels display significant correlations between IRAC 
colours, from 3.6 to S/xm. The third region includes the bot- 
tom left corner and shows IRAC colours vs. optical colours: 
correlations in this region are very weak if present at all. In 
star forming galaxies an exception to this is the [u] colour, 
which appears to (weakly) correlate with IRAC colours as 
well. This visual impression is quantitatively confirmed by 



the Spearman rank correlation coefficients, which we com- 
puted for each colour-colour distribution and are listed in 
the lower-right corner of each corresponding panel. 

From this analysis we can conclude that, on the local 
scales probed by our actual pixellation (from 200 pc up to 
< Ikpc when the effects of the smoothing are maximal), the 
shapes of the SED in a given galaxy in the optical range 
can be described, to first order, as a 1-parameter family 
and the same holds true for the IR range. The two parts 
of the SED shortwards and bluewards of 1.65/im appear to 
be largely uncorrelated. The optical SED is mainly driven 
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by light- weighted mean stehar age (and secondarily by dust 
extinction and metallicity, see e.g. Fig. Al in ZCR09 and 
MacArthur et al.|[200 ?) while the mid-IR SED probed by 
IRAC is determined by PAH features and hot dust, which 
trace on-going star formation (e.g. Calzetti et al.|200'7 ). This 
tells us that on local scales light- weighted mean stellar ages 
and actual on-going star formation are largely uncorrelated, 
most likely because stellar aging and star formation are sub- 
ject to very different time scales, namely ~ 10^~^ and ^10^ 



years, respectively (e.g. [Bruzual A. Charlot||1993 
cutt][l998t 



Kenni- 



It is also very notable that in star-forming galaxies the 
colours of the two shortest- wavelength IRAC bands at 3.6 
and 4.5/im correlate extremely well with those of the longer- 
wavelength bands, which are dominated by PAH features, 
but are very weakly correlated with the z and i bands which 
are very good tracers of the stellar mass (e.g. ZCR09). This 
is in agreement with the finding of Mentuch et al. (2010), 



who show, using pixel- by-pixel analysis, that the colour ex- 
cess at 3.6 and 4.5/im with respect to pure stellar emission 
extrapolated from J-band (1.25/im) is related to star for- 
mation. Therefore caution should be used when adopting 
images at 3.6//m alone to trace the stellar mass distribution 
in galaxies with no correction for the 'dust' contribution. 

3.2 Principal Component Analysis of resolved 
SEDs 

In this section we use principal component analysis (PCA) of 
the SEDs of all pixels, galaxy by galaxy, to show i) that the 
variation of the optical part and of the IR part of the SED 
are uncorrelated and ii) that these independent variations 
are spatially related to the structure of a galaxy. 

We refer the reader to Jolliffe ( 2002 ) for a thorough de- 



scription of principal component analysis. For the present 
purpose it is sufficient to say that the goal of PCA is to "re- 
duce the dimensionality of a data set consisting of a large 
number of interrelated variables, while retaining as much 
as possible of the variation present in the data set. This 
is achieved by transforming to a new set of variables, the 
principal components (PCs), which are uncorrelated, and 
which are ordered so that the first few retain most of the 
variation present in all of the original variables" ( Jolliffe| 
|2002|. In our specific application, for each pixel i of the 



iVpix pixels in a galaxy we consider the 9-dimensional SED 
vector given by the iif-band normalized logvfu^ i.e. SED^ = 

{\ogl>xfvx,i - ^OgiyHfuH,i}x=u,g,r,i,z,3.6A-5,5-7,8fj.m- The SCt 

of these A^pix SED vectors is then decomposed into principal 
components using the pcomp function in IDL. Each SED^ can 
be then represented as SED^ = SEDmean + X)j=i 9 cnj^^Ji 
where SEDmean is the arithmetic mean of SED^ taken over 
all pixels, PCj are the principal component vectors, which 
are obtained as eigenvectors of the correlation matrix of the 
dataset, and dijj are the corresponding coefficients for ith 
pixel and jth component. In Fig. [3] and ^ we illustrate the 
results of the PCA on our reference galaxies (corresponding 
figures for the complete sample are available in the online 
only appendix, figures CI to C5). Figure a) (top left) shows 



the pixel distribution of SEDs, with the arithmetic mean 
marked with blue crosses for reference. Figures b), c) and d) 
(to the right) display the pixel distribution in the SED space 
when only PCI, PC2 and PC3, respectively, are used in ad- 



dition to the mean, such that SED^,j = SEDmean + aijPCj, 
for j = 1, 2, 3 respectively. As apparent from these plots, the 
majority of the scatter is accounted for by the first two PCs. 
This holds true for the two examples shown here and in gen- 
eral for the full sample. The variance in each PC is given by 
the number in parenthesis next to the PC label in each plot. 
We find that typically ^90% of the variance is contained 
in the first 2 PCs (with PCI typically contributing at least 
roughly twice as much as PC2) and only 2 to 7% in PC3. 
The remaining variance is contributed in decreasing order by 
the other PCs, and is basically at the level of noise. While 
in general PCs are not expected to have any obvious physi- 
cal interpretation (this is certainly the case for PCi, i > 3), 
these plots show that PCI and PC2 mainly affect the IR and 
the optical part of the SED very differently, but in a very 
consistent way for different morphological types and levels 
of star formation activity. In particular, PCI accounts for 
most of the variations in the IR part of the SED and, in 
some cases, for part of the variations in the optical, whereas 
PC2 mainly reflects variations in the optical regime. This is 
seen more quantitatively with the Spearman rank correla- 
tion coefficients between colors and PC coefficients reported 
for each band at the top of figures b), c) and d), with these 
numbers indicating the relative contribution and direction 
to each band. This relation of PCI and PC2 holds true for 
five of the seven galaxies in our sampl^ with the only excep- 
tions being the inclined spiral NGC 3521 and the elliptical 
galaxy NGC 4552, which we discuss below in Sec. [3T7] The 
third principal component, though apparently dominated by 
variance in the optical bands for most galaxies, is only a low 
level contribution and is different enough between galaxies 
that it is of very uncertain physical interpretation and thus 
will not be discussed further in this work. 

The PCA demonstrates that the optical and IR parts 
of the SED are to large extent independent and their linear 
combination can reproduce pixels SED to high accuracy. As 
already noted in the previous section, this shows that on 
local scales (0.2 to < 1 kpc) the PAH emission that drives 
IR colours is not related to the combination of mean light- 
weighted age of stars and extinction, which is traced by the 
optical colours. Apparently any combination of PAH emis- 
sion and optical colours is allowed within a galaxy. 

The decomposition into PCs provides a very compact 
representation of the SEDs, which is essentially given by the 
two coefficients describing the intensity of the two main PCs. 
In Fig.|3]s^ and|4]s^ we show how the intensities of PCI and 
PC 2 affect different regions of the galaxies. The intensity 
of the red tint represents the intensity of PCI, while the 
blue one is for PC2. Regions dominated by either PCs in 
fact extend over a broad range of scales and trace a struc- 
ture in which spiral arms, bright star-forming knots therein, 
interarm regions, bulge can be identified. 

PCI (mainly IR-related) is typically strong (positive, 
enhanced IR emission) where bright star forming regions 
are found and is faint (negative, diminished IR emission) 



2 For NGC 4536 the role of the first two PCs is switched: PCI 
correlates best with optical colours, whereas PC2 antz-correlates 
with IR colours. However, modulo index exchange and sign of 
the correlation, the same observations and conclusions as for the 
other four "regular" galaxies apply to NGC 4536 as well. 
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Figure 3. SED pixel distribution and PCA for NGC4254. a) displays the distribution of all pixels in terms of if-band normalized SED. 
At each wavelength, corresponding to the ten optical-IR broad bands, the different colour intensities represent the number of pixels in 
bins of 0.02 dex of log ufiy, as shown by the key on the right side. The mean SED (i.e. the arithmetic mean of log ufiy for all pixels at 
each wavelength) is marked with blue crosses, b), c), d) are analogous to figure aj, except that the SED of each pixel is not the observed 
one, but the one reconstructed by using only the mean plus one single PC (1, 2 and 3 respectively) times the corresponding coefficient. 
The colour key is the same as in figure a). For each wavelength the Spearman rank correlation coefficient between the corresponding 
ffux and the strength of the PC is written at the top of the graph. This number indicates both the relative strength and sign of the 
PC contribution to that wavelength. The number in parenthesis beside the PC label is the fractional amount of variance contributed 
by the PC. Blue crosses are the mean SED, as in figure a). In e) we simultaneously map the intensity of PCI and PC2 on the galaxy 
image. In the RGB rendition, the intensity of PCI is mapped to the red channel (spanning the range corresponding to the range between 
the 5th and 95th percentile of the distribution in PCI coefficient) and PC2 is similarly mapped to the blue channel. The green channel 
intensity is just the mean of the blue and the red ones. Note that the range of PC coefficients is centered on and extends both positive 
(bright on the map) and negative (dark on the map). The overlaid contours correspond to different isophotal levels of Ho; emission: 
the five levels, ordered by increasing line thickness, correspond to logSB(HQ;) of —17.00, —16.25, —15.50, —14.75, —14.00 in units of 
erg s"-*^ cm~^ arcsec"^. 
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Figure 4. Same as Fig. [3] but for NGC4450. 
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in the bulge and in the galaxy outskirts without ongoing 
star formation. The very close link between star formation 
and PCI is analyzed in more detail in the following Sec. 
|3.3| and is expected from the observed relations between 
PAH, mainly responsible for enhanced emission in the IRAC 
bands, and SFR (see, e.g., jCalzetti et al. 2007). By over- 
laying Ha isophotal contours on the maps of Fig. [sjs^ and 
we show that PCl-bright regions are actually coincident 
with Ha-bright regions and therefore PCI can be broadly 
interpreted as an indicator of SFR density. On the other 
hand, PC2 (mainly optical-related) is strong (i.e. positive, 
enhanced optical luminosity with respect to NIR) where 
young, unattenuated stellar associations are found along the 
spiral arms, but also in the disk outskirts, as an effect of gra- 
dients in age (and possibly also in metallicity) , but is weak 



(i.e. negative, suppressed optical emission with respect to 
NIR) in the bulge, in the interarm space and in heavily at- 
tenuated regions in the arm^ This complex interplay be- 
tween regions of different physical conditions results in the 
observed lack of correlation between the optical and the IR 
part of the SED. As can be inferred from figures |3]s^ and 
|4]s^, this result appears to be robust against the choice of 
resolution scale, as long as a sufficient number of elements 
(pixels) can be analyzed and the limit at which the internal 
structure of the star forming regions start to be resolved is 
not reached (in fact, these scales are never reached in the 



^ Note that the interarm regions have a slightly negative PC2 
and a largely negative PCI, hence the resulting dark blue tint. 
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galaxies explored here as they are almost a factor of 10 below 
our resolution). 

It is notable that in the galaxies of later type (hence 
having higher specific SFR) PCI and PC2 appear to con- 
tribute very similar amount of variance in the [u] colour, 
which has, in fact, a very similar degree of correlation with 
the two main PCs. This can be easily accommodated in the 
interpretation given in the previous paragraph by consider- 
ing that the flux in u band is contributed both by star form- 
ing regions, hence correlated with PCI, and by the regular 
diffuse stellar population in relation to its age and relative 
attenuation, hence correlated with PC2. 



3.2.1 PC A of the "anomalous" 
NGC3521 



laxies NGC 4552 and 



Contrary to all other galaxies, the PCA of the elliptical 
NGC 4552 results in PCI being responsible for most (89%) 
of the total SED variance. All pixels in NGC 4552 have SEDs 
very close to a pure stellar continuum, with no evidence for 
dust/PAH emission and the IR part being consistent with 
a Ray leigh- Jeans tail. PCI affects the SED in the sense of 
changing the "curvature", i.e. of enhancing or decreasing at 
the same time all bands relative to H. PCI has a clear gra- 
dient, such that the central regions have redder optical- if 
colours and bluer H-IR colours, while the outskirts display 
the opposite. This kind of colour gradients are well known in 
the optical and NIR (see, e.g., | Michard|2000[ [2005|) and are 
usually well explained by metallicity gradients fKobayashi 
Arimoto|1999 i Henry & Worthey 1999), plus possibly minor 
contributions of age gradients as well. Our analysis shows 
that there is only one parameter driving the SED variations 
in NGC 4552, and this is likely to be linked to the metallicity 
of the stellar population. 

The case of NGC 3521 is opposite and much more com- 
plicated. First of all, the scatter in SED from pixel to pixel 
is evidently larger than in any other galaxy. The two main 
PCs of this galaxy contribute very similar amount of vari- 
ance (54 and 40% respectively), at odds with the fact that 
PCI is clearly the dominant contributor to the variance in 
all other galaxies. PCI appears to be more closely corre- 
lated to optical colours than PC2, which, as opposed, has 
stronger correlations with the IR bands; however PCI and 
PC2 do not separate optical and IR in a clear way. The spa- 



tial distribution of PCI and PC2 intensity (Fig. CI) shows 
a surprising result. In the inner parts PCI is approximately 
axissymmetric and reflects in part the spiral structure of the 
galaxy; beyond ^ 50 arcsec, the symmetry is broken, with 
the western part of the galaxy displaying more intense, posi- 
tive PCI, and the eastern part having mainly negative PCI. 
By looking at the colour-composite giH image of this galaxy 
(see ZCR09, Fig. 5), we see that this reflects the differential 
attenuation and the scattering of light in the inclined disc. 
Thick dust lanes on the near side (to the west) strongly at- 
tenuate all optical wavelengths, and even attenuate the H 
band as well, meaning that our normalization will be bi- 
ased for all these pixels. The opposite scenario is observed 
on the far side of the galaxy (to the east) where the dust 
back-scatters the light coming from the central regions, thus 
resulting in much bluer colours than expected from the in- 
trinsic stellar population. 

PC2 is clearly strong in the central concentrated bulge and 



in the very faint outskirts: this can be explained by the fact 
that PC2 anti-correlates with the colours which are most 
sensitive to current star-formation and young stars. 

In NGC 3521 the interplay between star-formation, stel- 
lar age and dust appears much more complex than in other 
galaxies and this would explain why the PCA in this case 
produces very different results than in other galaxies. As 
we show in the following sections, this is also reflected in 
NGC 3521 being an outlier in a number of other relations. 

3.3 Dependence of SEDs on surface brightness 
and star formation rate density 

The PCA performed in the previous section allows us to con- 
dense the properties of pixel SEDs into just two parameters 
(the intensity of PCI and PC2) instead of 9, thus making 
it much easier to study the dependence of such properties 
on other local parameters. In particular we analyze here the 
dependence of SED properties on stellar density and star 
formation rate density. 

To this goal, for each galaxy we bin all pixels in the 
2-dimensional space of log SB(if)-log SB(Ha) and compute 
the median value of the coefficients of PCI and PC2 in each 
bin. This is shown in the left and right diagrams, respec- 
tively, of Fig.[5]for NGC 4254 and Fig.|6]for NGC 4450 (fig- 
ures for the other galaxies are provided in the online ap- 
pendix [d] except for NGC 4552, which has no Ha imaging). 
Importantly, these quantities have direct physical interpre- 
tations, with Ha SB often used as a proxy for star formation 
rate density (Kennicutt 1998), while the if-band SB can 



be considered a rough proxy for stellar mass density (e.g. 
ZCR09). 

Regular face-on, late type, star forming galaxies, such 
as NGC 4254, display clear trends for PCI and PC2 intensity 
in this SB space. At any given SB(i7) PCI increases with in- 
creasing SB (Ha), while a reversed trend is seen with SB(i7) 
at given SB (Ha). The diagonal lines drawn in each figure 
mark the loci of constant SB(Ha)/SB(if). Constant PCI 
intensity levels are almost parallel to these lines, thus indi- 
cating that PCI is effectively a measure of SB(Ha)/SB(i7) 
or, in more physical terms, of the local specific star forma- 
tion rate sSFR = SFR/M^ ^ SB(Ha)/SB(H). On the other 
hand, this is not surprising, given the very tight correlation 
that we show to exist between PCI and the 8/im color (which 
is proportional to SB(8/im)/SB(if)) and the fact that the 
8/im emission is a very good tracer of SFR, similar to Ha. 
The PC2 intensity trends are somewhat more noisy, with 
PC2 decreasing at increasing SB(i7), while no significant 
trend with SB (Ha) is observed. Following the physical inter- 
pretation given in the previous section of PC2 as a tracer of 
optical colors and hence of mean stellar age (and dust extinc- 
tion) and considering that SB (if) is with good approxima- 
tion a monotonically decreasing function of radius, this PC2 
trend with SB (if) can be broadly interpreted as reflecting 
stellar age gradients in late-type galaxies (e.g. MacArthur 
eraL][2004| 



Although not as clear as for NGC 4254, PCI trends 
which can be interpreted as tracing sSFR are seen or hinted 
at also for earlier type galaxies, such as NGC 4450. Con- 
cerning PC2 it can be seen that three of the five "regular" 
galaxies (NGC 4254, 4321, 4450) display a trend of PC2 de- 
creasing at increasing SB (if), thus consistent with the hy- 
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Figure 5. Dependence of PCI and PC2 intensity on if-band and Ho; surface brightness in NGC4254. Left: median intensity of PCI 
for the pixels in bins of log SB(i^)-log SB(Hq;) (width is 0.075 and 0.15 dex respectively), colour coded according to the key on the right 
side. The diagonal lines are lines of constant SB(HQ;)/SB(iif) ^ sSFR. Bins with less than three contributing pixels are not displayed. 
Right: same as the left, but for PC2. 
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Figure 6. Same as Fig. [s] but for NGC4450. 



pothesis that PC2 reflects mean age gradients. However for 
the other two such a trend is not detected, possibly due to 
the presence of thick dust lanes (NGC4536) or very weak 
stellar population gradients (NGC4579). 



These trends of PCs with SB are not seen in NGC 3521 
and NGC 4552, whose first two PCs, for different reasons, 
do not have a straightforward physical interpretation, as al- 
ready discussed in Sec. |3.2.1| 



4 FROM INTERNAL TO GLOBAL SED 
CORRELATIONS 

In this section we show how to reconcile the lack of correla- 
tion between variations in optical and IR colours observed 
on local scales with the existence of such correlations be- 
tween integrated colours of galaxies. In Fig. [7| we combine 
the pixels from all seven galaxies in the sample and show the 
combined colour-colour plots, like those analyzed in Sec. [3] 
The coloured symbols identify the colours from the integral 
photometry of each galaxy, as indicated in the legend. Sym- 
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Figure 7. The distribution of pixels of all galaxies in the "colour-colour" spaces defined by the fluxes log ufi, normalized to if-band, 
in u, g, r, i, z, 3.6, 4.5, 5.7 and 8/im bands. The intensity of grey linearly traces the density of pixels at each point in the colour-colour 
spaces. The Spearman rank correlation coeflicient among the pixels in each panel is written at its bottom right corner. The histograms 
on the diagonal show the distribution of pixels in normalized log ufi, for each band as labeled on the y-axis. Coloured fllled dots indicate 
the integral colours of individual galaxies, according to the legend at the top. Colours are chosen such that the sequence from red to 
purple is sorted by the morphological type of the galaxy, i.e. from early- to late-type. 



bol colours from red to purple are used to represent "increas- 
ing" morphological types, from early to late. As for individ- 
ual galaxies, the optical colours of pixels are typically cor- 
related to very high degree, as are the IR colours. However, 
contrary to what we observe in individual galaxies, when the 
pixels from all galaxies are combined, the optical colours cor- 
relate with IR colours as well, albeit with a larger scatter 



than that observed between optical colours or between IR 
colours alone. The correlation is quantified in each panel by 
the Spearman rank correlation coefficient, which is typically 
> +0.8 for optical-optical and IR-IR colours, and > +0.5 for 
optical-IR colours. The global correlations in the optical vs. 
IR colours result from the fact that, although each galaxy 
covers a limited region in the colour-colour diagrams with 
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highly scattered clouds of points, the barycentres of the dis- 
tributions follow well defined relations. This can be seen in 
the colour-colour relations of the integrated points, where 
all galaxies, with the only notable exception of NGC3521 
(the blue point), follow very well defined sequences. In other 
words, we see that the shapes of the integrated SEDs are 
very correlated between galaxies: IR colours are very tightly 
linked to each other; by enhancing the IR flux relative to H, 
bluer optical- iif colours result, which, in turn, are very well 
correlated with each other. This is a natural consequence of 
galaxies' SFHs being typically smooth, such that a higher 
current rate of star formation (which causes enhanced IR 
flux) is generally accompanied by a higher abundance of 
young blue stars, formed in the recent past and already lib- 
erated from their dusty birth cloud ( [Chariot Fan||2000 ). 
The only outlier in our sample to most of these correlations 
is NGC3521, whose aspect and geometrical distribution of 
dust result in a complex and peculiar interplay between the 
attenuation and scattering of photons and hence in the ob- 



served anomalous colours (see also Sec. 3.2.1 and 5.3) 



We note that not only is the relation between IR colours 
extremely tight for the integrated SEDs and for the average 
pixel distribution, but also almost all pixels, no matter which 
galaxy they belong to, appear to follow extremely tight uni- 
versal relations, whose observed scatter is consistent with 
or just slightly larger than pure measurement errors. We 
analyze this in more detail in Sec. [5] 

The degree of correlation between the [S/xm] colour and 
other colours reaches a minimum with the SDSS z and i 
bands (Spearman rank correlation coeflficient of +0.31 and 
+0.51 respectively), with the correlation with the z band 
being the lowest for all colour combinations. This can be 
understood by considering that the [z] colour {[z] — [H]) 
and the [i] colour ([i] — [H]) are almost insensitive to stellar 
age, marginally sensitive to dust extinction and are mainly 
driven by metallicity (e.g. ZCR09), while the [8/im] colour 
is essentially a sSFR indicator (as we showed above), whose 
relation with metallicity is not obvious (see, e.g., [Calzetti 



et al. 2007 Gordon et al. 2008). Colours at shorter wave- 



lengths are increasingly related to stellar age and therefore 
it is not surprising that a better correlation with the [S/xm] 
colour is observed. This reaches a maximum with the cor- 
relation of the [8/im] colour with the [u] colour (for which 
a signiflcant correlation was already observed within some 
individual galaxies), as this colour (effectively u — H) can 
be also used as a rough proxy of sSFR. 

One last point to make about Fig. [7| is that in these 
colour-colour diagrams galaxies appear to be ordered in 
their integrated SEDs according to their morphological type 
(coded by the colours of the symbols). While this is not al- 
ways exactly the case when one looks at optical colours, the 
order is exact in the IR colours. [Li et al. ( 2007 ) also observed 
similar trends of S/xm (corrected for stellar contamination) 
vs. 3.6/im with morphology. This might be surprising since 
morphologies are determined at optical wavelengths, not IR. 
In physical terms, this implies that the morphology of galax- 
ies reflects primarily their sSFR, rather than the properties 
of their optical SED. 



5 MID-IR COLOUR RELATIONS 

5.1 A universal sequence of mid-IR colours and 
its predictive power 

In this section we analyze in more detail the relation between 
IR colours and demonstrate that the SEDs in this spectral 
window are a one-parameter family to very high accuracy. In 
the mid panels of Figures [sja^, h) and we plot the flux at 
3.6, 4.5 and 5.7 /xm, respectively, relative to the flux at 8 fim. 
as a function of the [8/im] — [H] colour. Greyscale levels rep- 
resent the log of pixel densities in the colour-colour spaces 
for the flve "regular" galaxies: NGC4552 and NGC3521 are 
analyzed separately, because of their anomalous natures (see 



also Sec. 5.3). The orange lines show the running percentiles 
for this distribution: 2.5, 16, 50, 84 and 97.2%, from the 
bottom to the top, respectively. The pixels of the elliptical 
NGC 4552 and the inclined dusty spiral NGC 3521 are repre- 
sented separately by the coloured points, with NGC 4552 the 
compact cloud of red points top left in each flgure (very low 
8 /im flux relative to H and to all other IRAC bands) , while 
NGC 3521 is the extended sequence of blue points partly 
offset from the sequence deflned by the main sample. These 
plots suggest that the ratio between 8 /xm and if-band lumi- 
nosities can be used to predict the luminosities in the other 
three IRAC bands with high accuracy. 

The bottom panels of Figures [sja^ to show the half- 
width of the 16-84% percentile range, a robust proxy for 
the r.m.s., of the pixels belonging only to the five "regular" 
galaxies. The horizontal dotted line corresponds to the max- 
imum r.m.s. that is expected for null intrinsic scatter and the 
minimum S/N of 20 in each band. At low [8/im] — [if] < —0.6, 
the scatter in the [3.6/im] — [8/im] and [4.5/im] — [8/im] colours 
is lower than the maximum possible contribution from mea- 
surement error, namely 0.02 dex or 5%, and the relation is 
close to linear. At higher [8/im] — [H] the relations fiatten and 
the scatter increases up to 0.035-0.040 dex (8-10%), thus re- 
quiring significant intrinsic scatter to be explained. As far as 
scatter is concerned, for the [5.7/im] — [8/im] colour the situa- 
tion is reversed: [8/im] — [H] can be used to predict this colour 
to typical accuracy better than 7% at [8/im] — [H] > — 1, but 
its performance is significantly worse (10-12% accuracy) at 
lower relative fiux intensity. The slope of the relation also 
changes significantly from almost —1 at low [8/im] — [H] to 
at high values. 

We fit fourth degree polynomials to the relations ob- 
served for the five "regular" galaxies. The coeflficients are 
given m Table [3] The fits are shown as solid black lines in 
the mid panels of the Figures [sja^ to c^, where they almost 
perfectly overlap to the running median lines. These rela- 
tions enable the generation of maps of the galaxies at 3.6, 
4.5 and 5.7 /im just from the i7-band and 8 /im images. An 
example of this is shown in Fig. |9] for NGC 4254; the images 
of the remaining six galaxies are provided for comparison in 
the online appendix (figures |E1| to |E6| ) . The three rows of 
images are for the IRAC 3.6, 4.5, and 5.7 /im bands respec- 
tively. The first image in each row is the original image as 
observed, the image in the middle is the one "synthesized" 
from the iif-band and 8 /im images using the fit relations 
given in Table [3] and the last image is the ratio between the 
synthetic image and the observed one. This figure, together 
with those relative to the other four "regular" galaxies in 



the online appendix (Fig. E2 E3 E4 E6), illustrates three 
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Figure 8. IRAC IR colour correlations, a), b) and c): mid panels show the distributions of pixels in the plane given by the flux at 3.6, 
4.5 and 5.7 /im (a, b and c respectively) relative to the flux at 8 /im vs. the flux at 8 /im relative to the flux in if-band (all fluxes are 
expressed in units of logufjy). NGC3521 and NGC4552 are represented by the blue and the red points respectively, while the density 
distribution of pixels belonging to the other flve galaxies are represented by the grey scale shading (log scaling). For reference the global 
colours of the seven galaxies are plotted as fllled circles, coloured according to the key in Fig. [T] The flve orange lines show the running 
median and 2.5, 16, 84 and 97.5 percentiles of the pixel distribution for the flve galaxies represented by the grey scale shading. The solid 
black curve is a 4th-degree polynomial fltting to these pixels, whose coeflicients are given in Table [3] The solid black straight line to 
the lower-left shows a slope of —1 to guide the eye. It is immediately clear that NGC3521 and NGC4552 are both outliers with respect 
to the rest of the sample. The points with error bars represent the pure stellar colours derived in Sec. |5.2| Theoretical predictions from 
CB07 for stellar population of ages > 2Gyr and ~ IGyr are marked with a x and a +, respectively, while the Starburst99 predictions 
adopted by Helou et al. (2004) are marked with solid horizontal lines. In the bottom panels we plot the r.m.s. (the half-widths of the 
16-84 percentile range) of the pixel distributions along the y-axis for the flve galaxies, as a function of [8/im] — [H]. The horizontal dashed 
line indicates the maximum r.m.s. possible due to measurement and photometric errors. In the top panels the local slope of the fltting 
polynomial is plotted. For reference, the slopes of (corresponding to constant flux ratio to 8/im) and —1 (corresponding to constant 
flux ratio to H) are marked as dotted lines, d): for each pixel we deflne the quantity "max Alogz^/zy" as the maximum of the absolute 
deviation between \ogufu observed and predicted based on the respective fltted polynomial functions, at 3.6)Lim, 4.5)Lim and 5.7 fim (see 
text for details). The fraction of pixels for which such a maximum deviation. A, is less than some difference, is plotted as solid lines, 
with blue showing the result for all galaxies, while black shows the result when the outliers NGC 3521 and NGC 4552 are excluded. The 
dotted lines show the fraction of if-band flux contained in pixels for which the maximum log ufu deviation is less than with blue for 
all and black excluding the outliers as for the other curves. 
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Figure 9. The images of galaxies (NGC4254 in this case) at 3.6, 4.5 and 5.7 fim can be reconstructed with great accuracy by combining 
the if-band and the 8/iin images only by means of the fitting functions provided in Tab. [3] On each of the three rows, for 3.6, 4.5 and 5.7 
fim from top to bottom, we reproduce, from left to right, the original image of the galaxy, the synthetic image reconstructed from if-band 
and the 8/im images, and the ratio between the synthetic and the observed images. While most of the pixels display typical deviations of 
a few per cent, consistent with or just slightly larger than photometric uncertainties, larger deviations up to a few 10% display a regular 
structure that hints at systematic departures from the universal fitting functions related to different local physical conditions. 



points: i) the overall accuracy of the synthetic images is on 
average very good, typically within a few per cent; ii) the 
same relations can be applied to all galaxies as no significant 
bias is seen from galaxy to galaxy; Hi) significant differences 
between the observed and the synthetic images are observed 
in some pixels, up to the level of 20-30%. The spatial dis- 
tribution of these pixels with the largest deviations is not 
random, hence not simply a consequence of noise, but rather 



reflects to some extent the structure visible in the intensity 
images. For instance, some low- SB regions appear to be bi- 
ased low in the synthetic 5.7 /xm image; arm and interarm 
regions appear to be differently biased with respect to each 
other in all three reconstructed bands. 

To better quantify point i) in the previous paragraph, 
we measure both the fractional area and H-hand luminos- 
ity occupied by the pixels for which the predicted IRAC 
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Table 3. Polynomial fit coefficients to the IR colour-colour rela- 
tions presented in Fig. [s] The y color named in column 1 can be 
obtained from [8/im] — [if] =xa.sy = cq-\-cix-\-C2X^ -\-c^x^ -{-c/^x^^ . 
Note that colours relative to [if], i.e. [3.6/im] — [ii], [4.5/xm] — [ii], 
[5. 7/im] — [ii], can be obtained from analogous polynomial approx- 
imations as a function of [8/i,m] — [H] by simply replacing ci with 
ci + 1. 



Colour 




CO 


ci 


C2 


C3 


C4 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


[3.6/im] 


- [8/im] 


-0.564 


-0.639 


0.420 


0.291 


0.079 


[4.5/im] 


- [8/im] 


-0.840 


-0.572 


0.471 


0.321 


0.091 


[5.7/im] 


- [8/im] 


-0.312 


-0.114 


0.201 


0.092 


0.060 



fluxes deviate from the observed ones by more than a given 
amount. For each pixel we determine the maximum abso- 
lute difference "max Alogz^/i^" between the predicted and 
observed pixel fluxes logu fi^ in the three bands, i.e. 3.6, 4.5, 
and 5.7/im. We plot the fraction of pixels for which "max 
A log z^/i/" is less than a given amount x and the fraction of 
luminosity in those pixels as a function of x in Fig. [s] c?^ as 
solid and dotted lines, respectively. Black lines are for the 
five "regular" galaxies only, while for the blue lines all seven 
galaxies are used. If we restrict the analysis to the "regular" 
galaxies only, in 80% of the pixels the maximum deviation 
between predicted and measured fluxes is < 0.052 dex, i.e. 
< 12.7%; 90% of pixels are in agreement within 0.069 dex, 
i.e. 17.2%. In terms of luminosity fraction, the measured ac- 
curacy is slightly worse, thus indicating that the main rea- 
son of disagreement between the observed colours and those 
predicted from our simple fitting formulae does not originate 
from noise only (or mainly): in fact this plot shows that, on 
average, the deviating pixels have if -band surface brightness 
higher than average. This confirms the point Hi) of the pre- 
vious paragraph, extracted from the visual analysis of Fig. 
[9] (and analogues in the online Appendix). 



5.2 Relative dust and star contributions and 
constraints on the purely stellar SED 

The IR colour-colour relations analyzed in the previous sub- 
section indicate that the knowledge of the flux ratio between 
the band mostly affected by dust and PAHs (8 /im) and the 
one including stellar emission only (if) is sufficient to deter- 
mine the full SED in between. To first order, we can repre- 
sent the IR SED as a simple linear combination of a stellar 
and a dust component, basically templates with constant 
contributions to each band: at high values of [S/nm] — [if] 
the IR SED is strongly dominated by the dust component, 
while the stellar continuum dominates at low values. The 
flux at the three intermediate wavelengths is expected to 
asymptotically approach a constant ratio to [8/im] going to- 
ward high [8/im] — [if], whereas a constant ratio to [if] is 
expected to be approached at low values. These two regimes 
translate into slopes of and —1 respectively in the colour- 
colour relations of Fig. [S] In particular, the slope of —1 
follows from the fact that for each [X] (3.6, 4.5 and 5.7 /xm) 
[X] - [8/im] = {[X]- [if]) - ([8/im] - [H]); as we approach 
the purely stellar regime, [X] — [if] should become asymp- 
totically constant and therefore [X] — [8/im] vs. [8/im] — [if] 



should be asymptotically approximated by a line of slope 
-1. 

It is therefore interesting to show the actual slopes as a 
function of [8/xm] — [if], which we plot in the upper panels 
of Fig. [s] a) to c). As far as 3.6 and 4.5 /xm are concerned, 
while slopes close to —1 exist for a wide range, with slopes 
a < —0.8 measured up to [8/im] — [if] = —0.5, a slope of 
a = — 1 is reached only at the very lowest [8/im] — [if] ra- 
tios. Past [8/im] — [if] ^ 0.0, a quickly approaches for 
both bands, though it never reaches in the actual data 
range for either. What these slopes imply is that, while the 
dust contribution is weak for most values of [8/im] — [if], the 
fluxes at 3.6 and 4.5 /im can be considered to be "pure" stel- 
lar continuum only for the very lowest [8/im] — [if] ratios; 
as soon as some PAH emission appears at 8 /xm, also 3.6 
and 4.5 /xm are affected. At [8/im] — [if] > —0.5 typical of 
star-forming galaxies (see Fig.[7|, these bands appear to be 
heavily contaminated by the PAH feature at 3.3 /xm, in the 
3.6 /xm band, and by hot dust and extended wings of neigh- 
bouring PAH features, in the 4.5 /im band (e.g.|Draine et al 



2007 Draine & Li 2007). This issue of dust contamination 



in the first two IRAC bands is also consistent with work by 
Meidt et al. (2011, ApJ submitted) and the increase in con- 
tamination with star formation was also found by[Mentuch| 
aI] ( [20T0| . 

The slope of the [5.7/im] — [8/im] vs. [8/im] — [if] relation 
indicates that the 5.7 fim follows the 8 fim emission much 
more closely (as expected). At the lowest [S/nm] — [if] ratios, 
slope does not even reach —1, indicating that nowhere in 
the five "regular" galaxies does the emission at 5.7 /xm arise 
purely from stellar continuum. The slope a quickly increases 
toward 0, although this value is only reached at the very 
highest [8/im] — [if] values: only there does the PAH/dust 
component fully dominate at 5.7 /xm, while at lower val- 
ues the contamination by stellar continuum is still signifi- 
cant. Interestingly, the scatter in the [5.7/im] — [8/xm] at the 
highest [8/im] — [if] values most likely reveals the maximal 
variation between the 6.6/im and 7.7/jm PAH features (and 
underlying continuum) that dominate these bands, at least 
in NGC 4254 which dominates these colours. This scatter is 
small, both relative to the expected noise (bottom panel of 
Fig. and to the variation observed in these PAH fea- 
tures between galaxies, as shown for the SINGS galaxies in 



Smith et al^ ( |2007[ especially Fig. 13). 

The exact knowledge of the purely stellar contribution 
to the galaxies' mid-IR SED is of great importance in order 
to gain a more accurate description of dust emission and 
its physics. As we showed in this section, this task is com- 
plicated by the fact that stellar and dust emission in this 
spectral region are always mixed to some degree. The com- 
monly adopted approach to decouple these two components 
heavily relies on models (both for stellar and for dust emis- 
sion 



or/and SED fitting techniques (see e.g. Helou et al 
2004 |da Cunha et al.|2008 ). It has to be stressed, however, 



that stellar population synthesis models in the mid-IR are 
still very uncertain and poorly constrained by observations. 
The diagnostic plots of Fig.jsja^, b) and c) allow us to derive 
empirical constraints on the purely stellar SED in the fol- 
lowing way. In the simple hypothesis that i) the SED at each 
pixel is just a linear combination of a constant stellar com- 
ponent and a constant dust component and ii) H band is not 
contaminated by dust emission, pure stellar emission would 
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Table 4. Pure stellar emission flux ratios derived from IR colour-colour relations presented in Fig. [s] A comparison to the predictions 
from stellar population synthesis models is given in the three bottom rows (see text for details). 



Slope from relation 


[3.6] - [H] 


[4.5] - [H] 


[5.7] - [H] 


[8] - [H] 


[4.5] - [3.6] 


[5.7] - [3.6] 


[8] - [3.6] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


[3.6/xm] - [8/^m] 


-0.742 


-1.044 


-1.194 


-1.491 


-0.301 


-0.453 


-0.749 


[4.5/xm] — [Syum] 


-0.740 


-1.046 


-1.178 


-1.384 


-0.306 


-0.437 


-0.644 


[5.7/j,m] — [8/j,m] 


-0.741 


-1.038 


-1.199 


-1.572 


-0.297 


-0.458 


-0.832 


mean 


-0.741 ± 0.001 


-1.043 ± 0.003 


-1.190 ± 0.009 


-1.482 ± 0.077 


-0.301 ± 0.004 


-0.449 ± 0.009 


-0.742 ± 0.077 


Helou et al. (2004) 










-0.225 


-0.399 


-0.635 


CB07 (> 2Gyr) 


-0.745 


-1.047 


-1.321 


-1.687 


-0.302 


-0.576 


-0.942 


CB07 (w IGyr) 


-0.644 


-0.893 


-1.158 


-1.514 


-0.249 


-0.514 


-0.870 



correspond, in principle, to the minimum possible value of 
[8/im] — [H] ([8/im] is chosen because it is the band most sen- 
sitive to dust emission and because [8/im] — [H] offers the 
largest wavelength leverage). However, empirically it is very 
diflacult to determine such a minimum because it is a limit 
which is never actually reached (no galaxy is fully dust-free) 
and because photometric errors and intrinsic SED variations 
of the stellar component (which are predicted by SPS mod- 
els) increase the uncertainties. We use the simple analytical 
formulae that fit the empirical relations of Fig.|8]a^, h) and 
to extrapolate the flux ratios to the limit of pure stellar 
emission, which we determine as the point ([S/xm] — [H]) at 
which the derivative of the relation is a = —1, according 
to the argument given previously. The determination of the 
purely stellar [8/im] — [H] depends on the relation one con- 
siders: for [3.6/im] — [8/im] and [5.7/im] — [8/im] the solution is 
actually found by extrapolation, whereas for [4.5/im] — [8/im] 
a = — 1 occurs within the range covered by the datcjj The 
variations in this "stellar point" among the three different 
determinations provides some direct empirical estimate of 
the uncertainty of this flux ratio. Flux ratios between the 
other bands are obtained by plugging the resulting "stellar" 
[8/im] — [H] colour into the corresponding 4th order poly- 
nomials from Table [3] The resulting three determinations of 
the stellar flux ratios are given in Table [4] along with their 
mean and standard deviation. Our determinations are ac- 
curate at < 0.01 dex (2%) level for all bands except [8/im], 
whose ratio with respect to shorter wavelength cannot be 
estimated better than 20%. We also report the IR flux ra- 
tios from SPS models: the wide-spread [Helou et al." (2004) 
values, which were computed using Starburst99 (Leitherer 
T999 



et al, 



and those we derive from a preliminary 2007 
version (CB07 hereafter) of the [Bruzual L Chariot] ( |2003 ) 
SPS models, which includes an improved treatment of TP- 
AGB stars. These models predict little dependence on stel- 
lar metallicity and age for ages larger than 2 Gyr, while 
variation up to 30% or more can be expected for younger 
populations. For this reason we report the ratios for simple 
stellar populations (SSPs) older than 2 Gyr and for an SSP 1 
Gyr old, both using solar metallicity and Salpeter IMF. The 
[Helou et al.| ([2004]) values typically over-estimate by ^ 20% 
the IR flux ratios to [3.6/im] with respect to our empirical 
estimates. By contrast, the old CB07 SSPs are in excellent 



^ This is possibly due to the fact that the 4.5/im band is the only 
band not directly including PAH features and therefore might 
reach a constant ratio to H at higher [8/im] — [H] than for the 
two other bands. 



agreement with observed estimates between H and 4.5/im, 
while it appears that fluxes at longer wavelengths are signif- 
icantly under-estimated by up to ~ 50%. In consideration of 
the apparent disagreement between SPS models, the empir- 
ical ratios derived from our analysis are of great importance 
to improve the SED modeling of this spectral region. 

5.3 Mid-IR colour relations of the "anomalous" 
galaxies NGC 4552 and NGC 3521 

The two galaxies that stand out due to their anomalous 
principal component analyses also have evident anomalies 
when we look at their IR colour-colour relations (i.e. red 
and blue clouds in Fig. [8|. 

As expected from the lack of current star-formation, 
the elliptical galaxy NGC 4552 has the lowest ratio of 8 /im 
luminosity to H and also the lowest ratios of luminosities 
in the three shorter IRAC bands with respect to [8//m], 
as one can see from the red points in Fig. [S] However, a 
substantial fraction of those points appear to lie above the 
extrapolated colour relations for all other galaxies. Stellar 
population synthesis models of these spectral regions are 
unfortunately still too inaccurate to provide an explanation 
for these deviations of the order of 20-30%. The observed 
phenomenon can be summarized by saying that the lumi- 
nosity at shorter IR wavelengths is (increasingly) suppressed 
with respect to longer wavelengths. As we verify using the 
stellar population synthesis models with dust attenuation of 
ZCR09, these mid-IR colours could be explained by invok- 
ing a substantial amount of reddening by dust, which can 
in principle move the data points above and to the right of 
the relations of Fig. [S] However, our maps show that the 
outer regions of the galaxy are the most affected by this 
effect and that there is a suppression of NIR emission also 
with respect to optical wavelengths. In addition, NGC 4552 
is not detected at 160 /xm with Spitzer ( [Draine et al.|2007| ), 
thus implying a negligible amount of dust. All these facts 
argue against dust reddening as a cause of the anomalous 
IR colours. We also explore the possibility that a radial vari- 
ation of the stellar mass function originates such an effect: 
using CB07 SPS models we flnd that for stellar populations 
older than a couple of Gyr (as it is the case for NGC 4552), 
the IR colours differ by less than a few per cent going from 
a Salpeter to a Chabrier IMF. According to the spectral 
indices measured by Trager et al. (2000) in different aper- 



tures and the Sauron index maps of Kuntschner et al. ( 2006 ), 



NGC 4552 is iron-enriched and a-enhanced to supersolar val- 
ues (by comparison with, e.g., the [Thomas et al.|2003| mod- 
els). [a/Fe] appears to moderately increase at large radii. 
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while the decrease in the strength of the iron indices toward 
larger radii hints at a moderate decrease in metallicity. Ac- 
cording to the CB07 SPS models, the decrease in metallicity 
could affect the mid-IR colours in the required direction. We 
note however that the anti-correlation between the strength 
of the iron indices and [a/Fe] ( [Thomas et al. 112003 ) makes 
this metallicity gradient quite uncertain. We suggest that 
the increase of [a/Fe] to values up to 0.35-0.40 dex higher 
than solar at large radii could be the reason of the anoma- 
lous (near)IR colours. A thorough investigation of the de- 
pendence of these colour on [a/Fe] and metallicity, which is 
beyond the scope of this paper, would be required to provide 
a robust explanation of the observed effects. 

As already discussed in Sec. |3. 2.1] anomalies in the SED 
of NGC3521 result from the complex effects of substantial 



El 



dust absorption and scattering. From Fig. |8p^ and Fig. 
(online appendix), there is a significant fraction of pixels in 
this galaxy in good agreement with the "universal" relations 
derived in the previous section. These pixels are located in 
the far (north-east) part of the inner disc, those apparently 
less affected by extinction and scatter. All other pixels ap- 
pear to form a sequence offset from the relations of Fig. [S] 
which underestimate the fluxes at 3.6, 4.5 and 5.7/im with 
respect to the observed ones. On the near (south-western) 
side of the disc the stellar light is heavily dust obscured, so 
much so that even the H-hand (and possibly the shortest 
IRAC bands) is attenuated, and this probably explains the 
offset of these pixels to the right (and up) of the relations. 
However, the reason for the pixels in the halo of the galaxy 
being similarly offset is unclear. It is possible that something 
similar to what happens in the outer regions of NGC 4552 
occurs here as well, as the continuity between the sequence 
of the red and the blue pixels of Fig. [S] suggests. 



6 SUMMARY AND CONCLUDING REMARKS 

Using mult i wavelength data (i^-band to 8/im) for seven 
nearby galaxies, resolved at scales of > 200 pc pix~^ with 
a Signal-to-Noise Ratio of > 10 — 20 for each pixel, we 
have demonstrated, both qualitatively and quantitatively, 
the disconnect that occurs between the optical and infrared 
colours relative to the near-IR if-band on the small scales 
within galaxies. A very important point in the analysis is 
the choice of the if-band point as the normalization of the 
pixel SEDs, as this band is the least affected by dust ab- 
sorption and emission at the same time and divides the por- 
tion of SEDs possibly affected by dust emission (longward 
of it) from the one dominated by stellar emission, possibly 
extincted by dust (shortward of it). Normalizing at other 
wavelengths where dust effects are more substantial would 
likely destroy the observed correlations by mixing different 
emission/absorption mechanisms. In fact, we have checked 
that using a slightly shorter wavelength for normalization, 
such as the z band, results in similar although more scat- 
tered correlations. On the other hand, if we move the nor- 
malization to longer wavelengths, we expect quite dramatic 
effects already at 3.6/im, where we have demonstrated the 
contamination by dust emission to be substantial. 

Qualitatively, the disconnect between the optical and 
infrared colours can be seen in colour-colour plots (fig- 
ures [l] [2j and Bl to B5) where correlations can be visi- 



bly seen between optical-optical colours, IR-IR colours, but 
not optical-IR colours. This holds true for all but two of the 
seven galaxies: NGC 4552, an elliptical galaxy whose colours 
at all wavelengths examined are dominated by stellar light 
from old stars, and thus all colours are correlated to a high 
degree; and NGC 3521, a more inclined late- type galaxy with 
significant scatter in its colours, with dust causing weird 
colour correlations due to the effects of attenuation on the 
near side of the galaxy, and scattering on the far side. 

The observed disconnect in colours can be investigated 
using Principal Component Analysis (PCA) of the distri- 
bution in if-band normalized Spectral Energy Distribution 
(SED) space of the pixels in each galaxy. This analysis picks 
out two distinct and uncorrelated components in each of the 
five "typical" galaxies; a component that dominates the vari- 
ance and is strongly correlated to the variations in the IR, 
and a component that is mainly correlated to the variations 
in the optical (see figures [3) [4] and [CTI 



to 



C5). These two 



components account for most (> 90%) of the observed vari- 
ance in colours across the galaxies, and show distinct spatial 
variations in galaxies, with the dominant primary compo- 
nent associated with the IR variance (PCI) being stronger 
in spiral arms and star forming regions, and the other op- 
tically dominated component (PC2) being enhanced in the 
interarm regions and outskirts of galaxies. 

These two principal components can be associated with 
physical properties within the galaxies by comparison with 
the surface brightnesses in Ha and in if-band, themselves 
proxies for the star formation rate (SFR) and stellar mass 
(M*). This comparison (figures [s] [g] and |Dl| to |D4[ ) reveals 
a strong correlation of PCI with increasing SB (Ha) /SB (if), 
an indicator of the specific star formation rate (sSFR), while 
PC2 variation correlates with SB (if), although with signif- 
icant scatter. 

Altogether, these correlations indicate that PCI is 
detecting regions where young stars are dominating the 
colours, and predominantly in the IR, while PC2 is picking 
up variations in the mean stellar age, which is associated 
with the stellar mass density, with the most massive regions 
generally the oldest on average. It is these two distinct com- 
ponents that both the colours and PCA technique reveal. 

Given the disconnect observed on the small scales, it 
is surprising that the global correlations of galaxy colours 
and types exist (e.g. Kennicutt 1998). However, as figure 
|7| reveals, while there are distinctions on the pixel scales 
(> 200pc), globally the mean SEDs and colours trace very 
well the type of the galaxy, with later-types being bluer and 
dustier (IR brighter) than early-types. Within each galaxy, 
the scatter (and disconnect of IR and optical) is around 
this mean SED. In more physical terms, we can interpret 
this as the result of sSFR and stellar mass density (stellar 
age) being correlated between galaxies, but not (necessarily) 
within galaxies. 

Interestingly, the relations between the mid-IR colours 
are even tighter when all seven galaxies in our sample are 
considered together. When an IR-IR colour (e.g. [3.6/im] — 
[8/im]) is plotted as a function of a IR-near-IR colour 
(i.e. [S/xm] — [H]) these tight relations are apparent (figure 
|8|. Simple 4th-order polynomial fits to these relations (re- 
sulting in the coefficients in Table [sj have a r.m.s. scatter 
around them which is, most of the time, less than the max- 
imum expected photometric uncertainties. These functional 



Break down of optical- IR relations 19 



fits are tight enougii tiiat tiie [8/im] — [H] colour images can 
be used as a predictive tool for the fluxes in the three other 
IRAC bands with deviation of only a few percent typically 
(figure [9] and figures 



El 



to 



E6) 



What these relationships reveal is the varying contribu- 
tion of dust to the IRAC bands. The if-band is purely stel- 
lar, with no dust contribution, while the S/xm IRAC band is 
dominated by the 7.7/im poly cyclic aromatic hydrocarbon 
(PAH) complex, with only a weak stellar contribution, and 
thus the [8/im] — [H] colour is a measure of relative dust 
contribution to the pixel SEDs. The colour-colour diagrams 
of figure [S] therefore show the contribution of dust and stars 
in each band, going from pure stellar (where the slope is 
— 1) to dust dominated (where the slope is 0). By extrapo- 
lating the function fits to where the slope becomes —1 we 
are able to empirically determine the pure stellar colours in 
the IRAC bands (table |3| , at least for the galaxies in our 
sample. The plots clearly show the need to correct for dust 
when using the IRAC bands as tracers of stellar mass (see 
e.g. Meidt et al. 2011 for such a correction). Interestingly, 
these colours match that of the 2007 models of !Bruzual M 
[Chariot (2003) (not calibrated for these wavelengths) for old 
(> 2Gyr) stars at the shortest wavelengths, but is much red- 
der (flatter SED) at the long (5.7 and S/xm) bands and bluer 
(steeper SED) than that suggested by |Helou et al. (2004) in 
these bands. Another interesting note is that the scatter of 
pixels in the stellar and dust regimes of these relations gives 
an indication of the intrinsic scatter in the stellar SED or the 
"dust" SED respectively on the scales explored here. This 
scatter is surprisingly low (figure [sja^ -6^, lower panels), given 
the known contribution of various stellar evolutionary stages 
to the if-band, though the scatter is somewhat larger for 
dust given the observed variation in PAH features and hot 



dust contribution in galaxies with Spitzer-IRS (e.g. [Smith 
et al.[[2007| 



In Summary, while on global scales the IR and optical 
colours of galaxies are correlated with their star formation 
history and galaxy- type, the correlation of IR and optical 
colours breaks down internally in galaxies on local scales. 
The IR dominated and optical dominated components can 
be spatially associated with regions and with physical pro- 
cesses, with the IR excess coming from regions of high spe- 
cific star formation rate, while the optical excess is broadly 
associated to low if-band surface brightness (stellar mass 
density) regions and to unobscured young stellar associa- 
tions along spiral arms. The disconnect of optical and IR 
colours is because of this spatial distinction of the uniform 
underlying older stellar population, and the more concen- 
trated, stochastic star forming regions. However, this discon- 
nect of optical and IR colours in galaxies is around a mean, 
and the mean colours of galaxies follow the well-observed 
trend of bluer optical and stronger IR in later type galaxies. 
The IRAC colours within galaxies follow each other closely, 
with a tight enough relationship that a IRAC-near-IR colour 
can be used to predict the fluxes in the other IRAC bands. 
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A//i,m 

Figure Al. Images in the eleven bands used for the analysis in this work for NGC3521. The greyscale indicates the log of the surface 
brightness in each band in MJySr"-*^ as defined by the colour bar to the right of each panel, except for the Ha image which is in log 
ergcm"^ s"-*^ arcsec"^. In the lower right panel we show the distribution of logufi, SEDs of the pixels for the galaxy, normalized to the 
if-band. The Ho; fiux is also represented in this plot (at A = 0.6563/im) using a fiux density corresponding to the fiux divided by a 
passband of lA. The colour scale indicates the number density of pixels (in bins of 0.02 dex in normalized \ogvfu) as labelled in the key 
to the right. The blue symbols indicate the arithmetic mean of the pixels ('x' symbol) and the luminosity-weighted mean of the pixels 
(i.e. integrated SED of the galaxy, symbol). 



APPENDIX A: IMAGES IN ELEVEN BANDS AND SEDS FOR THE FULL SAMPLE 
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Figure A3. As in Fig[AT] but for NGC4321. 
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Figure A4. As in Fig[Al] but for NGC4450. 
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Figure A5. As in Fig |Al| but for NGC4536. The stripe of missing data going through the nucleus of the galaxy is due to the detector 
saturation in the IRAC bands. 
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Figure Bl. As in Fig^ but for NGC3521. 



APPENDIX B: COLOUR CORRELATIONS FOR THE FIVE GALAXIES NOT IN THE MAIN PAPER 
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Figure B2. As in FigjT] but for NGC4321. 
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Figure B3. As in Fig^ but for NGC4536. 
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Figure B4. As in Fig^ but for NGC4552. 
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Figure B5. As in Fig[Tl but for NGC4579. 



Break down of optical-IR relations 33 




ARA/arcsec d) i 

A/ fim 

Figure CI. As in Fig[3] but for the galaxy NGC3521. 



APPENDIX C: PCA DECOMPOSITION FOR THE FIVE GALAXIES NOT IN THE MAIN PAPER 
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Figure C3. As in Fig[3] but for the galaxy NGC4536. The stripe of missing data going through the nucleus of the galaxy is due to the 
detector saturation in the IRAC bands. Note that in this case, contrary to the other four "regular" galaxies, PCI is mainly correlated to 
the optical colours and PC2 is anti-correlated with the IR colours. Modulo this exchange, main conclusions remain unaffected. 
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Figure C4. As in Fig[3] but for the galaxy 



NGC4552. 
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Figure C5. As in Fig[3] but for the galaxy NGC4579. 
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Figure Dl. As in Figjs] but for NGC3521. 



PCI 




-0.5 0.0 0.5 1.0 1.5 
Log SB(H)/MJy Sr"^ 



CP 



X 

m 

CP 

o 




-0.5 0.0 0.5 1.0 1.5 

Log SB(H)/MJy Sr"^ 



Figure D2. As in Figfs] but for NGC4321. 



APPENDIX D: DEPENDENCE OF PRINCIPAL COMPONENTS ON SURFACE BRIGHTNESS AND 
STAR FORMATION FOR THE FOUR GALAXIES NOT IN THE MAIN PAPER 
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Figure D3. As in Figjs] but for NGC4536. Note that in this case, contrary to the other four "regular" galaxies, PCI is mainly correlated 
to the optical colours and PC2 is anti-correlated with the IR colours. Modulo this exchange, main conclusions remain unaffected. 
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Figure D4. As in Fig[5] but for NGC4579. 
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Figure El. As in Fig[9] but for NGC3521. 
APPENDIX E: SYNTHETIC IRAC IMAGES FOR THE SIX GALAXIES NOT IN THE MAIN PAPER 
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Figure E2. As in Fig[9] but for NGC4321. 
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Figure E3. As in Fig[9] but for NGC4450. 
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Figure E4. As in Fig[9] but for NGC4536. 
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Figure E5. As in Fig[9] but for NGC4552. 
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Figure E6. As in Fig[9] but for NGC4579. 



